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Abstract 

A model of dark matter and dark energy based on the concept of gravitational polarization is 
investigated. We propose an action in standard general relativity for describing, at some effective or 
phenomenological level, the dynamics of a dipolar medium, i.e. one endowed with a dipole moment 
vector, and polarizable in a gravitational field. Using first-order cosmological perturbations, we 
show that the dipolar fluid is undistinguishable from standard dark energy (a cosmological constant 
CN | A) plus standard dark matter (a pressureless perfect fluid), and therefore benefits from the successes 

of the A-CDM (A-cold dark matter) scenario at cosmological scales. Invoking an argument of "weak 
clusterisation" of the mass distribution of dipole moments, we find that the dipolar dark matter 
reproduces the phenomenology of the modified Newtonian dynamics (MOND) at galactic scales. 
The dipolar medium action naturally contains a cosmological constant, and we show that if the 
model is to come from some fundamental underlying physics, the cosmological constant A should be 
of the order of a^/c 4 , where ao denotes the MOND constant acceleration scale, in good agreement 
with observations. 
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I. INTRODUCTION 

In the current concordance model of cosmology (the A-CDM scenario, see e.g. [l|) based 
on Einstein's general relativity (GR), the mass-energy content of the Universe is made of 
roughly 4% of baryons, 23% of cold dark matter (CDM) and 73% of dark energy in the form 
of a cosmological constant A. The dark matter accounts for the well-known discrepancy 
between the mass of a typical cluster of galaxies as deduced from its luminosity, and the 
Newtonian dynamical mass [2J . The model has so far been very successful in reproducing the 
observed cosmic microwave background (CMB) spectrum [3j and explaining the distribution 
of baryonic matter from galaxy clusters scale up to cosmological scales by the non-linear 
growth of initial perturbations [4j. Although the exact nature of the hypothetical dark 
matter particle remains unknown, super-symmetric extensions of the standard model of 
particle physics predict well-motivated candidates (see [5[] for a review). Simulations suggest 
some universal dark matter density profile around galaxies [6|. However, in that respect, the 
CDM hypothesis has some difficulties [7], |8| at explaining in a natural way the distribution 
and properties of dark matter at galactic scales. 

The modified Newtonian dynamics (MOND) was proposed by Milgrom [9|, llC], bj| to ac- 



count for the basic phenomenology of dark matter in galactic halos, as evidenced by the flat 
rotation curves of galaxies, and the Tully-Fisher relation [12j between the observed luminos- 
ity and the asymptotic rotation velocity of spiral galaxies. However, if MOND serves very 
well for these purposes (and some others also [§]), we know that MOND does not fully ac- 
count for the inferred dark matter at the intermediate scale of clusters of galaxies [13] , |l4J, |l5 



In addition, MOND cannot be considered as a viable physical model, but only as an ad- 
hoc — though extremely useful — phenomenological "recipe". In the usual interpretation, 



MOND is viewed (see [16| for a review) as a modification of the fundamental law of gravity 
or the fundamental law of dynamics, without the need for dark matter. The relativistic 
extensions of MOND, of which the Tensor- Vector-Scalar (TeVeS) theory [14j, |l7|, ll8( is the 



prime example, share this view of modifying the gravity sector, by postulating some suple- 
mentary fields associated with the gravitational force, in addition to the metric tensor field 
of GR (see [19| for a review). Recently, such modified gravity theories have evolved toward 



Einstein-aether like theories XM 21, 22,^3, 24 



Each of these alternatives has proved to be very successful in complementary domains of 
validity: the cosmological scale (and cluster scale) for the CDM paradigm and the galactic 
scale for MOND. It is frustrating that two successful models seem to be fundamentally 
incompatible. In the present paper we shall propose a third approach, which has the potential 
of bringing together the main aspects of both A-CDM and MOND in a single relativistic 
model. Namely, we keep the standard law of gravity, i.e. GR and its Newtonian limit, but we 
add to the distribution of ordinary matter some specific non-standard form of dark matter 
(described by a relativistic action in usual GR) in such a way as to naturally explain the 
phenomenology of MOND at galactic scales. Furthermore, we prove that this form of dark 
matter leads to the same predictions as for the A-CDM cosmological scenario at large scales. 
In particular, we find that the relativistic action for this matter model naturally contains 
the dark energy in the form of a cosmological constant A. Thus, our model will benefit from 
both the successes of the A-CDM scenario, and the MOND phenomenology. 



The model will be based on the observation [25|, |26| that the phenomenology of MOND can 



be naturally interpreted by an effect of "gravitational polarization" of some dipolar medium 
constituting the dark matter. The effect can be essentially viewed (in a Newtonian-like 



interpretation [251 ]) as the gravitational analogue of the electric polarization of a dielectric 
material, whose atoms can be modelled by electric dipoles, in an applied electric field [27|. In 
the quasi-Newtonian model of |25j the gravitational polarization follows from a microscopic 
description of the dipole moments in analogy with electrostatics. It was shown that the 
gravitational dipole moments require the existence of some internal non-gravitational force 
to stabilize them in a gravitational field. Thanks to this internal force, an equilibrium 
state for the dipolar particle is possible, in which the dipole moment is aligned with the 
gravitational field and the medium is polarized. The MOND equation follows from that 
equilibrium configuration. However the model [25J cannot be considered as viable because it 
is non-relativistic, and involves negative gravitational-type masses (or gravitational charges) 
and consecutively a violation of the equivalence principle at a fundamental level. 

In a second model |26j we showed that it is possible to describe dipolar particles con- 
sistently with the equivalence principle by an action principle in standard GR. The action 
depends on the particle's position in space-time (as for an ordinary particle action) and 
also on a four-vector dipole moment carried by the particle. The particle's position and 
the dipole moment are considered to be two dynamical variables to be varied independently 
in the action. Furthermore, a force internal to the dipolar particle was introduced in the 
form of a scalar potential function (say V) in this action. The potential V depends on some 
adequately defined norm of the dipole moment vector. Because of that force, the particle is 
not a "test" particle and its motion in space-time is non-geodesic. The non-relativistic limit 



of the relativistic model 26] was found to be different from the quasi-Newtonian model [2 



(hence the two models are distinct) but it was possible under some hypothesis to recover 



the same equilibrium state yielding the MOND equation as in [25J. However the relativistic 
model [26j, if considered as a model for dark matter, has some drawbacks -- notably the 
mechanism of alignement of the dipole moment with the gravitational field is unclear (so 
the precise link with MOND is questionable), and the dynamics of the dipolar particles in 
the special case of spherical symmetry does not seem to be very physical. 

In the present paper, we shall propose a third model which will be based on an action 
similar to that of the relativistic model [26] but with some crucial modifications. First we 
shall add, with respect to [26J , an ordinary mass term in the action to represent the (inertial 
or passive gravitational) mass of the dipolar particles. Second, the main improvement we 
shall make is to assume that the internal force derives from a potential function in the action 
(call it W) which depends not on the dipole moment itself as in [26j but on the local density 
of dipole moments, i.e. the polarization field. In this new approach we are thus assuming 
that the motion of the dipolar particles is influenced by the density of the surrounding 
medium. This is analogous to the description of a plasma in electromagnetism in which the 
internal force, responsible for the plasma oscillations, depends on the density of the plasma 
(cf. the expression of the plasma frequency [27I]). 1 Because the action [given by (12. 2p with 
(12. 7p below] will now depend on the density of the medium, it becomes more advantageous 
to write it as a fluid action rather than as a particle action. 

This simple modification of the model, in which the potential W depends on the po- 
larization field, will have important consequences. First of all, the relation with the phe- 
nomenology of MOND will become clear and straightforward. Secondly, we shall find that 
the motion of dipolar particles in the central field of a spherical mass (in the non-relativistic 



1 In the quasi-Newtonian model [25] the dipolar medium was formulated as the gravitational analogue of a 
plasma, oscillating at its natural plasma frequency. 



limit) makes now sense physically. The drawbacks of the previous model [26(] are thus cured. 
Last but not least, we shall find that the model naturally involves a cosmological constant. 
Then, with the equations of motion and evolution (and stress-energy tensor) derived from 
the action, we show the following: 

1. The dipolar fluid is undistinguishable from standard dark energy (a cosmological con- 
stant) plus standard CDM (say a pressureless perfect fluid) at cosmological scales, i.e. 
at the level of first-order cosmological perturbations. 2 The model is thus consistent 
with the observations of the CMB fluctuations. However, the model should differ from 
A-CDM at the level of second-order cosmological perturbations. 

2. The MOND phenomenology of the flat rotation curves of galaxies and the Tully-Fisher 
relation is recovered at galactic scales (for a galaxy at low redshift) from the effect of 
gravitational polarization. There is a one-to-one correspondence between the MOND 
function (say /i = 1 + x) an d the potential function W introduced in the action. 

3. The minimum of the potential function W is a cosmological constant A. We find 
that if W is to be considered as "fundamental", i.e. coming from some fundamental 
underlying theory (presumably a quantum field theory), the cosmological constant 
should be numerically of the order of Oq/c 4 , where a^ denotes the MOND constant 
acceleration scale. 

A relation of the type A ~ ciq/c 4 between a cosmological observable A and a parameter ao 
measured from observations at galactic scales is quite remarkable and is in good agreement 
with observations. More precisely, if we define the natural acceleration scale associated with 
the cosmological constant, 

■* = 5F)/f (L1) 

then the current astrophysical measurements yield ao — 1.3 oa- The related numerical 
coincidence a ~ cH was pointed out very early on by Milgrom [9|, |l_0|, 111]. The near 



agreement between ao and oa has a natural explanation within our model, although the exact 
numerical coefficient between the two acceleration scales cannot be determined presently. 

Since the present model will not be connected to any (quantum) fundamental theory, it 
should be regarded merely as an "effective" or even "phenomenological" model. We shall 
even argue (though this remains open) that it may apply only at large scales, from the 
galactic scale up to cosmological scales, and not at smaller scales like in the Solar System. 
However, this model offers a nice unification between the dark energy in the form of A and 
the dark matter in the form of MOND (both effects of dark energy and dark matter occuring 
when gravity is weak). Furthermore, it reconciles in some sense the observations of dark 
matter on cosmological scales, where the evidence is for the standard CDM, and on galactic 
scales, which is the realm of MOND. It would be interesting to study the intermediate scale 
of clusters of galaxies and to see if the model is consistent with observations. Such a study 
should probably be performed using numerical methods. 



2 Note however that while in the standard scenario the CDM particle is, say, a well-motivated supersym- 
metric particle (perhaps to be discovered at the LHC in CERN), in our case the fundamental nature of 
the "dipolar particle" will remain unknown. 



The plan of this paper is as follows. In section [TT] we present the action principle for 
the dipolar medium, and we vary the action to obtain the equation of motion, the equation 
of evolution and the stress-energy tensor. In section IIHI we apply first-order cosmological 
perturbations (on a homogeneous and isotropic background) to prove that the dipolar fluid 
reproduces all the features of the standard dark matter paradigm at cosmological scales. We 
investigate the non-relativistic limit of the model in section HVl and show that, under some 
hypothesis, the polarization of the dipolar dark matter in the gravitational field of a galaxy 
results in an apparent modification of the law of gravity in agreement with the MOND 
paradigm. Section |V] summarizes and concludes the paper. The dynamics of the dipolar 
dark matter in the central gravitational field of a spherically symmetric mass distribution is 
investigated in appendix lAl 



II. DIPOLAR FLUID IN GENERAL RELATIVITY 

A. Action principle 

Our model will be based on a specific action functional for the dipolar fluid in standard 
GR. This fluid is described by the four-vector current density J M = an 11 , where u^ is the four- 
velocity of the fluid, normalized to g^u^u" = — 1, and where a = *J—g IMV J tJ 'J v represents 
its rest mass density 3 In this paper we shall conveniently rescale most of the variables used 



in [26j by a factor of 2m, where m is the mass parameter introduced in the action of [26 



Hence we have a = 2mn, where n is the number density of dipole moments in the notation 



of 26 . The above current vector is conserved in the sense that 



V = 0, (2.1) 

where V M denotes the covariant derivative associated with the metric g^. Our fundamental 
assumption is that the dipolar fluid is endowed with a dipole moment vector field £ M which 
will be considered as a dynamical variable. We have £ M = 7r M /2m where 7r M is the dipole 
moment variable used in [26j (hence £ M has the dimension of a length) . 

Adopting a fluid description of the dipolar matter rather than a particle formulation as 
in [26J, 4 we postulate that the dynamics of the dipolar fluid in a prescribed gravitational 
field g^ u is derived from an action of the type 



S= / d 4 a;v^L[^,e,e,^], (2-2) 




where g = det^j,), the integration being performed over the entire 4-dimensional manifold. 
This action is to be added to the Einstein-Hilbert action for gravity, and to the actions of 
all the other matter fields. The Lagrangian L depends on the current density J M , the dipole 



3 Greek indices take the space-time values /U, is, . . . = 0, 1, 2, 3 and Latin ones range on spatial values i, j, . 
1,2,3. The metric signature is (—,+,+,+). The convention for the Riemann curvature tensor R^^ 



the same as in [28]. Symmetrization of indices is {[iv) = \{pv + v[i) and (if) = ^(ij + ji). In sections HT1 
and IIIII we make use of geometrical units G = c = 1 . 
4 The fluid action is obtained from the particle one by the formal prescription ^2 J dr — * J d 4 x ^—gn, 

where the sum runs over all the particles, and n is the number density of the fluid. 



moment vector £ M , and its covariant derivative £ M with respect to the proper time r (such 
that dr = ^—g^dx^dx"), which is defined using a fluid formulation by 

? = -£- = u'V^" , (2.3) 

dr 

and where D/dr is denoted by an overdot. In addition, the Lagrangian depends explicitly on 
the metric g^ v which serves at lowering and raising indices, so that for instance £ M = gw^- 
We shall consider an action for the dipolar medium similar to the one proposed in [26], 
with however a crucial generalization in that the potential function therein, which is sup- 
posed to describe a non-gravitational force internal to the dipole moment, will be allowed 
to depend not only on the dipole moment variable £ M , but also on the rest mass density of 
the dipolar fluid a. More precisely, we shall assume that the potential function W in the 
action depends on the dipole moment £ M only through the polarization, namely the number 
density of dipole moments, that is defined by 

IF = a^ , (2.4) 

or equivalently IP = nir^ in the notation of [26]. The dynamics of dipolar particles will 
therefore be influenced by the local density of the medium, in analogy with the physics of 
a plasma in which the force responsible for the plasma oscillations depends on the density 
of the plasma J27|]. Our assumption is that W is a function solely of the norm II_l of the 
projection of the polarization field (12.41) perpendicular to the velocity, namely 



n_L = yjg^nini = ^imp . (2.5) 

Here, the orthogonal projection of the polarization vector reads 11^ = _L^ U u , with the 
associated projector defined by i.^ = g^ u + u^u v . Similarly, we can define £^ = i_£J £" and 
its norm £j_ so that the (scalar) polarization field reads 

II_l = a£ ± . (2.6) 

The chosen dependence of the internal potential on II j_ will result in important differences 
and improvements with respect to the model of [261 ] . 
Our proposal for the Lagrangian of the dipolar fluid is 



a 



■i-V(^-^)(«"-e) + ^e 



- W(U ± ) , (2.7) 



where the two dynamical fields are the conserved current vector J M = au^ and the dipole 
moment vector £ M . The fourth term is our fundamental potential which should in principle 
result from a more fundamental theory valid at some microscopic level. The third term in 
(j2.7p is the same as in the previous model 26j and clearly represents a kinetic-like term for 



the evolution of the dipole moment vector. This term will tell how this evolution should 
differ from parallel transport along the fluid lines. The second term in (12. 7p (also the same 
as in [26|]) is made of the norm of a space-like vector and is inspired by the known action 
for the dynamics of particles with spin moving in a background gravitational field [29J . The 
motivation for postulating this term is that a dipole moment can be seen as the "lever arm" 
of the spin considered as a classical angular momentum (see a discussion in [26j]). 

6 



Finally, we comment on the first term in (12. 7p which is a mass term in an ordinary sense. 
The dipolar fluid we are considering will not be purely dipolar (or mostly dipolar) as in the 
previous model [26] but will involve a monopolar contribution as well. Here we shall thus 
have some dark matter in the ordinary sense. The mass term in (12.71) has been included for 
cosmological considerations, so that we recover the ordinary dark matter component at large 
scales (see section UTT|) . However, one can argue that the presence of such mass term a is not 
fine-tuned. Indeed, this term corresponds to the simplest and most natural assumption that 
the relative contributions of this mass density and the second and third terms in (12.71) are 
comparable. In addition, we notice that a = 2m n corresponds to the inertial mass density 
of the dipole particles in the quasi- Newtonian model [25j , so it is natural by analogy with 
this model to include that mass contribution in the action. Notice however that, even if 
the dipolar fluid is endowed with a mass density in an ordinary sense, its dynamics is well- 
defined only when the dipole moment is non-zero. Indeed, we observe that the Lagrangian 
(12. 7p becomes ill-defined when £ M = since the second term in (12 .7p is imaginary. 



B. Equations of motion and evolution 

In order to obtain the equations governing the dynamics of the dipolar fluid, we vary the 
action (12.21) [with the explicit choice of the Lagrangian (12.71) ] with respect to the d yna mical 
variables £^ and J M . The calculation is very similar to the one performed in [26J, but 
because of the different notation adopted here for rescaled variables (e.g. £ M = 7r /i /2m), and 
especially because of the more general form of the potential function, we present all details 
of the derivation. Varying first with respect to the dipole moment variable £ M , the resulting 
Euler-Lagrange equation reads in general terms 5 

in which the partial derivatives of the Lagrangian in ( 12.21) are applied considering the four 
variables £ M , £ M , J^ and g^ v as independent. For the specific case of the Lagrangian (12.71) . 
we get what shall be interpreted as the equation of motion of the dipolar fluid in the form 

K* = -T» , (2.9) 

in which the left-hand-side (LHS) is the proper time derivative of the linear momentum 6 

K ti = e + k». (2.10) 



o 



We write the Euler-Lagrange equation in this particle-looking form to emphasize the fact that the action 
(|2.7p is a particle (or fluid) action. Of course, this equation is equivalent to the usual field equation 

„ / dL \ dL 



() 



\dV^J d^ 
The present notation is related to the one used in [26[ by if = P^ /2m, k^ = p^ /2m, JF^ = F^/m 



(and f = ■K ti /2m). The quantity called A in 26|] is now denoted S in order to avoid confusion with the 
cosmological constant. 



Here, we introduced like in 26| a special notation for a four-vector k p which is space-like, 



whose norm is normalized to k^k,, = 1, and which reads 



fc" 



■;//' 



with S 



i-2^e + ^e 



(2.H) 



The space-like four- vector k^ will not represent the linear momentum (per unit mass) of the 
particle — that role will be taken by K^ which, as we shall see, will normally be time-like, 
see (I2.20aj) below. The quantity H has an important status in the present formalism because 
it represents the second term in the Lagrangian ([2.7]) and we shall be able to set it to one 
in section III CI as a particular way of selecting some physically interesting solution. On the 
right-hand-side (RHS) of (12.91) . the force per unit mass acting on a dipolar fluid element is 
given by 

^ = n^W n± , (2.12) 

in which we denote the unit direction of the polarization vector by II j_ = IIj_/ITl — 6jV6_l 
and the ordinary derivative of the potential W by W n = dW/dllj.. The "internal" force 
(12. 12)) being proportional to the space-like four- vector £^ = _L^ £", we immediately get the 
constraint 

u lt J rfl = Q. (2.13) 

We now turn to the variation of the action with respect to the conserved current J M = au^ 
(hence we deduce a = \J—J V J V and u^ = J^/cr). The general form of the Lagrange equation 
for the conserved current density reads (see e.g. |30j) 7 



D f dL 



dr V d J* 



u v V, 



dL 
~dT v 



(2.14) 



For the case at hands of the Lagrangian ( 12.71) . we get the following equation, later to be 
interpreted as the evolution equation for the dipole moment, 



£> = -v M (w - n ± w n± ) - R^ puX u P CK x 



(2.15) 



A new type of linear momentum fi M - - having the same meaning as in [26|] - has been 
introduced and defined by 



o// = j' />■" with j' = a 1 ' { i + !&,r+£±w ni ) - n,i" p' vim) 



The Riemann curvature term in the RHS of (12. 15)) represents the analogue of the coupling 
to curvature in the Papapetrou equations of motion of particles with spin in an arbitrary 
background 3l|. The complete dynamics and evolution of the dipolar fluid is now encoded 
into the equations (12.91) and (12.151) . Such equations constitute the appropriate generalization 



7 This can alternatively be written with ordinary partial derivatives as 



<m 



dL 



for the case of a density-dependent potential W, and in fluid formulation, of similar results 



m 



26 



Notice that by contracting (I2.15P with J M , the second term in the RHS of (12.151) cancels 
because of the symmetries of the Riemann tensor, and we get 



D 

cb 7 



J ^ = — (>V-n ± W ni ). (2.17) 



One can readily check that this constraint (12.171) can alternatively be derived from the other 
equation (12.91) together with the definition of f2 M in (12.161) . On the other hand, contracting 
( 12. 9p with u^ yields u^K^ = 0, which according to the definition of K^ 1 , leads to the other 
constraint 

«^[(S-1)F]=0. (2.18) 

This constraint can be viewed as a differential equation for the variable 5. 



C. Particular solution of the equations 



Following 26J, we shall solve the constraint (12.181) with the most obvious and natural 
choice of solution that 

5 = 1. (2.19) 

We shall see that this choice greatly simplifies the other equations we have. In particular, 
we are going to prove that the equations of motion ( 12. 9j) and equations of evolution ( 12.151) . 
when reduced by the condition H = 1, finally depend only on the space-like component 
of the dipole moment that is orthogonal to the velocity, namely £j\ so that the time-like 
component along the velocity, i.e. u u £ u , will have no physically observable consequences 
(actually, in that case this unphysical component turns out to be complex [26j). 

The structure of the subsequent equations and the physical properties of the model will 
heavily rely on the condition 5 = 1. Note that we could regard this condition not as a choice 
of solution but rather as a choice of theory. Indeed, we are going to pick up the simplest 
theory out of a whole set of theories in which 5 could have some non trivial proper time 
evolution obeying (I2.18p . Actually, we can view the choice 5 = 1 as an elegant way to impose 
into the Lagrangian formalism the condition that in fine the only physical component of the 
dipole moment should be £j\ namely the one perpendicular to the four- velocity field. We 
can imagine that it would be possible to impose the same physical condition in a different 
way, for instance by using Lagrange multipliers into the initial action. For exemple, in 



TeVeS [14], [17|, [18| or in Einstein- aether gravity [20j, [21|, [22|, [23|, [24J, a dynamical time-like 
vector field whose norm is unity is introduced by this mean. However, the present situation 
is different because our final physical vector £^ is space-like. 

When the condition (12.191) holds, the two linear momenta (I2.10p and (12.161) simplify 
appreciably and we obtain 

K? = u" , (2.20a) 

W = v?(\ + e±W n± ) + ±»t± . (2.20b) 

We see that the linear momentum K^ is finally time-like. These expressions depend only on 
the orthogonal component £j\ and we denote £^ = D£^/dr. The equations of motion and 
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evolution take now the simple forms 



VT 



0" 



-n^w ni , 



1 



a 



vM w - n , w 



v pM 



&# 



pu\ 



Pn,* 



wu 



(2.21a) 
(2.21b) 



Finally, the whole dynamics of the dipolar fluid only depends on the space-like perpendicular 
projection £j* of the dipole moment. 



D. Expression of the stress-energy tensor 



We vary the action (j2.2p with respect to the metric g^ to obtain the stress-energy tensor. 
We must first consider the general case where 5 is unconstrained, and then only on the result 
make the restriction that 5=1. We properly take into account the metric contributions 
coming from the Christoffel symbols in the covariant time derivative £ M by using the Palatini 
formula 32] . We are also careful that while the dipole moment £ M should be kept fixed during 



the variation, the conserved current J M will vary because of the change in the volume element 
y/— gd A x. Instead of J 11 , the relevant metric-independent variable that has to be fixed is the 
"coordinate" current density defined by J p = sf—g J M . Straightforward calculations yield 
the expression of the stress-energy tensor for an action of the general type (j2.2p . We find 



dg^u V 9 J" J d¥ 

V di p d£ v ) d£ v) 



^2.22) 



in which we denote dL/d^ p = g px dL/d^ x . The partial derivatives of the Lagrangian are 
performed assuming that its "natural" arguments J M , £ M , £ M and g^ are independent. The 
application to the particular case of the Lagrangian (J2.7P gives, for the moment for a general 
value of 5, 



T i» = _ g »v f W _ n ± W n± ) + Qbf) _ y p ([^K^ - K p ^] J u) ) 



(2.23) 



In the second term of (J2.23P we see that the linear momentum f2 M is related to the monopolar 
contribution to the stress-energy tensor, while the other linear momentum K p parametrizes 
the dipolar contribution in the third term. Comparing with equation (2.14) of [26], we 
observe that a new term, proportional to the metric g^, has been introduced. This term 
will clearly be associated with a cosmological constant, and we shall discuss it in detail 
below. One can readily verify that the conservation law V V T^ U = holds as a consequence 
of the equation of conservation of matter (I2.ip . and the equations of motion and evolution 
(P and fl2~ToD . for general 5. 

In the next step we reduce the expression ( I2.23P by means of the condition 5 = 1 and get 



T^ = -W gv + <j[ u»u v + f ± W n J."" + u 



(p 



s c 



_L 



V, 



ftvP - u p i { l 



J"A . (2.24) 



Again we notice that this expression depends only on the perpendicular projection £^ of the 
dipole moment. 
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It will be useful in the following to decompose the stress-energy tensor (12.241) according 
to the general canonical form 

T^ = r u^u v + V IT + 2 Q^u v) + TT , (2.25) 

where r and V represent the energy density and pressure, where the "heat flow" Q p is 
orthogonal to the four-velocity, i.e. u^Q^ = 0, and the symmetric anisotropic stress tensor 
YT V is orthogonal to the four-velocity and traceless, i.e. uJD iV = and E^ = 0. We get 



r = u p u a T pa , (2.26a) 

1 
3 
Q lt = -±Jiu T'" i (2.26c) 



V=-± p(T T p<T , (2.26b) 



while the anisotropic stress tensor is obtained by subtraction. In the case H = 1 where the 
dipolar fluid is described by the stress-energy tensor ( 12.241) we find that the energy density, 
pressure, heat flow and anisotropic stress tensor read respectively 

r = W-n ± W n± + p, (2.27a) 

V = -W + ^U ± W Jl± , (2.27b) 

Q" = <rg + n ± w n± ^ - n x ± v x u p , (2.27c) 

TT = Q J."" - e ± C ± ) n ± W n± , (2.27d) 

where we denote £^ = £±/£±, and where we introduced for future use the convenient notation 

p = a - V A n^ . (2.28) 

By contrast to ordinary perfect fluids, the characteristic feature of the dipolar fluid is the 
existence of non- vanishing heat flow Q^ and anisotropic stresses Y^ v . Furthermore, we notice 
that the energy density r involves (via p) a dipolar contribution given by — V^II^. That 
contribution will play the crucial role, as we will see in section IIV[ when recovering the 
phenomenology of MOND. 

III. COSMOLOGICAL PERTURBATIONS AT LARGE SCALES 

We are going to show in this section that the model of dipolar dark matter [i.e. based on 
the action (12.21) and (12. 7p . with equations of motion reduced by the condition 5 = 1] contains 
the essential features of standard dark matter at cosmological scales. We shall indeed prove 
that, at first order in cosmological perturbations, it behaves like a pressureless perfect fluid. 
Furthermore, we shall see that the dipolar fluid naturally contains a cosmological constant 
(the interpretation of which will be discussed below), and is thus supported by the obser- 
vations of dark energy. The model is therefore consistent with cosmological observations of 
the CMB fluctuations. 
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A. Perturbation of the gravitational sector 

We apply the theory of first-order cosmological perturbations around a Friedman- 
Lemaitre-Robertson- Walker (FLRW) background. For every generic scalar field or com- 
ponent of a tensor field, say F, we shall write F = F + 5F, where the background part F is 
the value of F in a FLRW metric, while 5F is a first-order perturbation of this background 
value. 

The FLRW metric is written in the usual way in terms of the conformal time 77, such that 
dt = adr] where a(r]) is the scale factor and t the cosmic time, as 

ds 2 = g^ dar'W = a 2 [-drj 2 + 7y dxW] . (3.1) 

Here 7^ is the metric of maximally symmetric spatial hypersurfaces of constant curvature 
K = or K = ±1. The perturbed FLRW metric ds 2 = g^ v dx^dx u will be of the general 
form [33] 

ds 2 = a 2 [-(1 + 2A) dry 2 + 2hidrj dx l + (7^ + h^) dzW] . (3.2) 



Making use of the standard scalar- vector-tensor (SVT) decomposition [34|, |35(, the metric 
perturbations hi and h^ are decomposed according to 

hi = DiB + B t , (3.3a) 

hij = 2C 7ij + 2D i D J E + 2D {l E j} + 2E ij , (3.3b) 

where Di denotes the covariant derivative with respect to the spatial background metric 7^. 
The vectors B l and E % are divergenceless, and the tensor E 1 ^ is at once divergenceless and 
trace-free, i.e. 

DiB* = D t E l = , (3.4a) 

D { E ij = Ei = 0. (3.4b) 

Spatial indices are lowered and raised with 7^ and its inverse 7 jfc . From these definitions, 
one can construct the gauge-invariant perturbation variables 

$ = A + (B' + KB) - {E" + HE') , (3.5a) 

* = -C - U{B - E') , (3.5b) 

Qi^Ei-Bi, (3.5c) 

together with Eij which is already gauge- invariant. The prime stands for a derivative with 
respect to the conformal time rj, and Ti = a' /a denotes the conformal Hubble parameter. 
We shall also use the alternative definition for a gauge-invariant gravitational potential 

*^- c -(l)'=* + * + (0- (36) 

B. Kinematics of the dipolar fluid 

The four-velocity of the dipolar fluid is decomposed into a background part and a per- 
turbation, u^ = u^ + 5u^. We have both g fiu ll^u l/ = —1 and g^ u u^u v = — 1. The background 
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part is supposed to be co moving, that is u l = 0. This defines a zeroth order in the per- 
turbation. In a FLRW background this means that it will satisfy the background geodesic 
equation u = 0. With standard notations, we have 



^ = -(1,0), 
a 


(3.7a) 


5u» = -{-A,F), 
a v ' 


(3.7b) 


r ill be written as u M = u^ + 5m m , with 




u^ = a (-1,0), 


(3.8a) 


S UfM = a(-A,fii + hi). 


(3.8b) 



The velocities of all the other fluids (baryons, photons, neutrinos, . . . ) are decomposed in a 
similar way. The perturbation of the three- velocity j3 l is split into scalar and vector parts, 

P 1 = D l v + v i with Dp* = , (3.9) 

and we introduce the gauge-invariant variables describing the perturbed motion, 

V = v + E', (3.10a) 

Vi = Vi + Bi. (3.10b) 

The dipolar dark matter fluid differs from standard dark matter by the presence of the 
dipole moment £^ (satisfying u^^_ = 0) carried along the fluid trajectories. For the dipole 
moment we also write a decomposition into a background part plus a perturbation, namely 
£j_ = £j_ + 5£j\ However, because a non- vanishing background dipole moment would break 
the isotropy of space, and would therefore be incompatible with a FLRW metric, we must 
make the assumption that the dipole moment is zero in the background, so that it is purely 
perturbative. Hence, we pose 

S = 0, (3.11a) 

5^={0,X l ), (3.11b) 

where A* represents the first-order perturbation of the dipole moment. Beware of our nota- 
tion for which A* is a vector living in the background spatial metric 7^. Thus the covariant 
components of the dipole moment perturbation are #£j_„ = (0,a 2 Aj) where Aj = jijX^. Note 
that there is no time component in the dipole moment perturbation because of the constraint 
Un£j_ = which reduces to u^S^ = at linear order. Like for the three-velocity field j3 l in 
(13. 9p . we split A* into a scalar and a vector part, namely 

\ i = D i y + y i with Dtf = . (3.12) 

However, unlike for v and v\ we notice that y and y % are gauge- invariant perturbation 
variables. This is because the background quantity is zero, £j_ = 0, hence the perturbation 



8£± is gauge-invariant according to the Stewart- Walker lemma [36|, l37J . 
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C. Cosmological expansion of the fundamental potential 

The next step is to make more specific our fundamental potential function W(IIj_) enter- 
ing the Lagrangian (|2.7jl . Such function should be a "universal" function of the polarization 
of the dipolar medium, described by the polarization scalar field 

n ± = a^. (3.i3) 

Now, we have seen that in cosmology there is no background (FLRW) value for the dipole 
moment, hence the background value of the polarization field is zero: n_i_ = 0. In linear 
perturbations, the polarization is expected to stay around the background value. Therefore, 
it seems physically well-motivated that the value II j_ = corresponds to a minimum of the 
potential function W, so that n_i_ does not depart too much from this background value, at 
least in the linear perturbation regime. We therefore assume that W(ITjJ is given locally 8 
by an harmonic potential of the form 

w(n ± ) = w + ^w 2 ni + o (ny , (3.14) 

where Wo and W2 are two constant parameters, and we pose Wi = 0. For linear pertur- 
bations, because II j_ = 8H± is already perturbative, we shall be able to neglect the higher 
order terms 0(ITjJ in (13.141) because these will contribute to second order at least in the 
internal force (j2.12p . Inserting the ansatz (13.141) into (12.121) we obtain 

j^ = w 2 u^ + o{u 2 ± ). (3.15) 

We asserted in the previous section that the background motion of the dipolar fluid is 
geodesic, i.e. u — 0. This is now justified by the fact that the force H3. 151) drives the 
non-geodesic motion via the equation of motion (j2.21a|) . hence since this force vanishes in 
the background, the deviation from geodesic motion starts only at perturbation order. 

In the present model the coefficients Wo, W2, ... of the expansion of our fundamental 
potential W(II_|_) are free parameters, and therefore will have to be measured by cosmological 
or astronomical observations. First of all, it is clear from inspection of the action (12.71) . or 
from the general decomposition of the stress-energy tensor [see ( I2.27al) and (I2.27bj) ]. that 
Wo is nothing but a cosmological constant, and we find 

W = A (3.16) 

The coefficient Wo is thereby determined by cosmological measurements of "dark energy". 
As we shall show in section HV] the next two coefficients W2 and W3 will be fixed by requiring 
that our model reproduces the phenomenology of MOND at galactic scales [8[ , and we shall 
find that W2 = 47r and W3 = 327r 2 /ao where a® is the constant MOND acceleration scale. 

Hence, in this model the cosmological constant A appears as the minimum value of the 
potential function W, reached when the polarization field is exactly zero, that is on an 
exact FLRW background (see Fig. [3]). Thus, it is tempting to interpret A as a "vacuum 
polarization", i.e. the residual polarization which remains when the "classical" part of 



The domain of validity of this expansion will be made more precise in section IIVB1 
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the polarization II^ — ► 0. Of course our model is only classical, hence there is no notion 
of vacuum polarization which would be due to quantum fluctuations. However, we can 
imagine that the present model is an effective one, describing at some macroscopic level a 
more fundamental underlying quantum field theory (QFT) in which there is a non-vanishing 
vacuum expectation value (VEV) of a quantum polarization field giving rise to the observed 
cosmological constant |38j. Then, the constant Wo would play the role of the VEV of this 
hypothetical quantum polarization field in such a more fundamental QFT. 



D. Perturbation of the dipolar fluid equations 

As for the four-velocity w M = TZ M + <5w M , we consider a linear perturbation of the rest 
mass energy density of the dipolar fluid according to a = a + 5a. The conservation law 
V M (cr« /i ) = reduces in the case of the background to 

a' + 3Ha = 0, (3.17) 

hence a evolves like a~ 3 . Concerning the perturbation, we define a = a (1 + e) so that the 
rest mass density contrast reads 

5a 

e = — . (3.18) 

a 

This quantity is not gauge-invariant, and one can associate with it in the usual way a 

gauge-invariant variable by posing 

—I r~i 

e F = e-—— = e + 3C, (3.19) 

a n 

with the index F standing for "flat slicing" . Alternatively, it is possible to introduce other 
gauge-invariant variables, like for example 

£ N =e-3H(B-E') = e F + m, (3.20) 

where the index N stands for "Newtonian". For the linear perturbation, the conservation 
law V At (crM /i ) = gives the gauge-invariant equations 

ep + AV = 0, (3.21a) 

e' N + AV = 3V', (3.21b) 

where A = r jijD l D : ' denotes the usual Laplacian operator. In the following we shall choose 
to work only with the flat-slicing variable e-p. 

According to f)2.21al) . the motion of the dipolar fluid obeys the equation w M = — JP\ A 
straightforward calculation yields the gauge-invariant expression for the four-acceleration, 

u^ = — (o, £>*($ + v + nv) + v l ' + nv^ . (3.22) 

On the other hand, the force is given by (I3.15P at first-order in the perturbation, in which we 
can use Uj_ = (0, a\ l ) to this order, with A* = D l y + y 1 . Hence, in terms of gauge-invariant 
quantities, the scalar and vector parts of the equation of motion read 

V' + HV + § = -Anaa 2 y, (3.23a) 
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Vj + HVi = -47caa 2 yi . (3.23b) 

Here we are anticipating on the results of the section [IV] and have replaced the constant W2 
in the expression of the force (13.151) by its value Att determined from the comparison with 
MOND predictions. 

If there was no dipole moment (i.e. y = y l = 0), we would recover the standard geodesic 
equations for a perturbed pressureless perfect fluid (see e.g. [33j), and according to (13.23bl) . 
the vector modes would satisfy (aVi)' = 0, and therefore vanish like a~ l . In contrast with the 
standard perfect fluid case, the dipolar fluid may have non-vanishing vector modes because 
of the driving term proportional to yi. Equation (I3.23ap clearly shows that the scalar modes 
are also affected by a non-zero dipole moment. 

The equation of evolution of the dipole moment was given by (12.21bl) . Now, f2 M reduces 



to £j_ + v? at first perturbation order, hence the evolution equation gives at that order 

g + if = -ftE^iW , (3.24) 

where R is the Riemann tensor of the FLRW background. By easy calculations we find 
for the derivatives of the dipole moment variable 

Q = -(o,\ i ' + H)A, (3.25a) 

£> = 1 (0, A'" + n\ 1 ' + H'X 1 ) . (3.25b) 

The scalar and vector parts of the equation of evolution are thus given by 

j/' + ny' = -(V' + HV + ^), (3.26a) 

y'S + Hy'^-iVj + HVi). (3.26b) 

Notice that the equation for the vector modes can be integrated, giving the simple relation 

yi + Vi = ^, (3.27) 

where s, is an integration constant three-vector. 

A comment is in order at this stage. Recall that we have included in the original La- 
grangian (12.71) a mass term in the ordinary sense, with the most natural value of the mass 
density simply given by o. This choice was made having in mind the physical analogy with 
the quasi-Newtonian model [25j where a = 2m n represented the inertial mass of the dipolar 
particles. Now we can see on a more technical level that such mass term is in fact essential 
for the model to work properly. If this mass term was set to zero in the action, then the 
RHS of both equations (13.26aj) and (13.26bl) would be zero. We would then find that y' and 
y\ vanish like a -1 , so that the dipole moment would in fact rapidly disappear or at least 
become non-dynamical, and the whole model would turn out to be meaningless. 

Combining the equations of motion (I3.23P and the evolution equations (13.261) . we obtain 
some differential equations for the scalar and vector contributions y and y l of the dipole 
moment A* = D l y + y l , which turn out to be decoupled from the equations giving V and V\ 
and to be exactly the same, viz 

y" + Hy' -4naa 2 y = 0, (3.28a) 
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y" + Hy' l -4naa 2 y i = 0. (3.28b) 

We find it remarkable that the dipole moment decouples from the other perturbation vari- 
ables so that its evolution depends in fine only on background quantities, namely a and the 
scale factor a. Since the equations for the scalar and vector modes are the same, we have 
also the same equation for the dipole moment itself, 

X^ + n\' i -4:7iaa 2 \ i = 0. (3.29) 

Clearly, the solutions of (13.291) behave typically as increasing and decreasing exponentials 
moderated by a cosmologial damping term TC \[. We can also write this equation in terms 
of the cosmic time t = J adrj, namely 9 

\ i + 2H\ l -4na\ i = 0, (3.30) 

where H = a/a = a' is the usual Hubble parameter. We find that the equations (I3.29P or 
( I3.30p are the same as the equation governing the growth of the density contrast of a perfect 
fluid with vanishing pressure for the sub-Hubble modes (say k 3> H) and when we neglect 
the contribution of other fluids; see (I3.50p below. In particular this means that like for the 
case of the density of a perfect fluid there is no problem of divergence (i.e. blowing up) of 
the components of the dipole moment A, between, say, the end of the inflationary era and 
the recombination. We can thus apply the theory of first-order cosmological perturbations 
even for the components of the dipole moment itself, which should stay perturbative. 

Notice that the value of the coefficient W2 = 47r used in (I3.29P or (13.301) . which makes 
such equations identical with the equation of growth of cosmological structures in the stan- 
dard CDM scenario, will only be determined in section HVl from a comparison with MOND 
predictions. There is thus an interesting interplay between the cosmology at large scales 
and the gravitational physics of smaller scales. 10 



E. The perturbed stress-energy tensor 

Consider next the stress-energy tensor of the dipolar fluid, that we decomposed as (I2.25J) 
with the expressions (I2.27l) - (l2.28p . At first perturbation order, these expressions reduce to 

r = W + p, (3.31a) 

V = -Wo , (3.31b) 

Q" = -fo,W\ i 'Y (3.31c) 

W u = 0, (3.31d) 

together with 

p = a(l + e-D i X i ). (3.32) 



9 In this equation, the dot stands for a derivative with respect to the coordinate time t, and not the proper 

time t as everywhere else. 
10 Actually the coefficient Air in (|3.29|) could be changed if we had assumed a mass term in the action (|2.7p 
different from a (say 2a or cx/2). The simplest choice we have made (for different reasons) that a is the 
correct mass term in the action corresponds also to the usual-looking evolution equation (|3.29[) . 
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We first note that part of the dipolar medium is actually made of a fluid of "dark energy" 
satisfying p^ e = — -Pde = VVo = A/87T where A is the cosmological constant. Accordingly, we 
shall write the decomposition 

T^ = T^ + TZ, (3.33) 

where the stress-energy tensor associated with the cosmological constant is denoted by T^, 
and where the other part represents specifically a fluid of "dark matter" whose stress-energy 
tensor is T^. Their explicit expressions read 

?T = ->Vo<r, (3.34a) 

TZ = p u»u u + 2 Q^u u) . (3.34b) 

Note that the dark matter part of the dipolar fluid, which may be called dipolar dark matter, 
has no pressure P, no anisotropic stresses £ M ", but a heat flow Q M given by f!3.31c|) and an 
energy density p given by (13.321) . or alternatively 

p = a{l+e -Ay). (3.35) 

The background energy density is simply given by the background rest mass energy density, 
p = a, and the corresponding energy density contrast is 

5=^r = £- Ay. (3.36) 

It differs from the rest mass energy density contrast e because of the internal dipolar energy. 
Like for e, one can construct several gauge-invariant perturbations associated with 5. We 
shall limit ourselves to the flat-slicing (F) one defined by (recall that y is gauge- invariant) 

5 F = S + 3C = e F - Ay, (3.37) 

and whose evolution equation is 

S' F + AV + Ay' = . (3.38) 

Similar gauge- invariant density contrast variables are also defined for the other fluids. Next, 
we split the dark matter stress-energy tensor f!3.34bl) into a background part plus a linear 
perturbation, namely T^ = T dm + ST^, and find 

TZ, = pu< l u u , (3.39a) 

STZ = 5pu^u u + 2p5u^u" ) + 2Q^u u) . (3.39b) 

We made use of the fact that the heat flow Q^ is already perturbative to replace the four- 
velocity in the last term by its background value. 

We are now going to show that the dipolar dark matter stress-energy tensor is undistin- 
guishable, at linear perturbation order, from that of a perfect fluid with vanishing pressure. 
To this end, we introduce the effective perturbed four-velocity 

bW = Su" + ^. (3.40) 

P 
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5^ = - 



Dipolar dark 
matter worldline 
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5*? l = 



Effective perfect 
fluid worldline 



FIG. 1: Sketch of the equivalence at first order of cosmological perturbations between dipolar dark matter 
and an effective perfect fluid. The dipolar dark matter has a four- velocity u^ = u^ + Su^, and follows a 
non-geodesic motion driven by the internal force T^, namely ii^ = —T^ L . One can construct from u^ and 
the heat flux Q^ an effective four-velocity u M = vr + Su^ satisfying a geodesic motion, i.e. u = 0. This is 
the four-velocity field of the effective perfect fluid associated with dipolar dark matter. 



Notice that u h 



u^+5u^ is still an admissible velocity field because Su° = —A/ a by virtue of 
the transversality property w M Q M = 0. The perturbed part of the dark matter stress-energy 
tensor (13.39bj) can then be written in the form 



0J dm 



5pu»u u + 2-p5u^u u) 



(3.41) 



which, together with (j3.39a|) . is precisely the stress-energy tensor of a perfect fluid with 
vanishing pressure P, vanishing anisotropic stresses S Miy , and a four-velocity field u^ = 
TZ M + 5u^. Using the definition ( 13.401) of the perturbed four- velocity Su* 1 , with the explicit 
expression of the heat flow fl3.31cj) . one can check that this perfect fluid consistently follows 
a geodesic motion, i.e. Su = 0. 

More explicitly, we can write the latter effective perturbation of the four-velocity in the 
standard form 8u^ = a~ 1 {—A 1 (3 l ) 1 and find that the effective ordinary velocity reads 



P i = P i + X 



(3.42) 



which can be viewed as a modification of the space-like component of the dipolar dark matter 
four-velocity This allows one to build a new four-velocity which would be tangent to the 
worldline of the effective perfect fluid (cf. Fig. [1]). In terms of scalar and vector parts, if we 

write Pi = D^v + v^ then 



v + y , 

vt + y'i ■ 



(3.43a) 
(3.43b) 
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Like for the perturbed four- velocity <5w M , we can introduce the gauge- invariant variables 

V = v + E' = V + y', (3.44a) 

Vi^Vi + B^Vi + y'i. (3.44b) 

In terms of the gauge-invariant variables V, Vi and Sp, the dipolar dark matter fluid equations 
(13T23J) and ([338]) finally read 

V' + HV + $ = Q, (3.45a) 

Vl + HVi = 0, (3.45b) 

5' F + AV = . (3.45c) 

These are precisely the standard evolution equations of a perfect fluid with no pressure and 
no anisotropic stresses (see e.g. [331] ). 

To summarize, we have proved that at first order of perturbation theory — and only at 
that order — the dipolar fluid behaves exactly as ordinary dark energy (i.e. a cosmological 
constant) plus ordinary dark matter (i.e. a perfect fluid). If we specify the background rest 
mass energy density a so that f2d m = 87i1jo/3Hq ~ 0.23 today as evidenced in cosmological 
observations, we can assert that the first-order cosmological perturbation theory with the 
dipolar fluid described by the stress-energy tensor (13.33H — (13.3411 will lead to the same pre- 
dictions than the standard A-CDM scenario — and is therefore consistent with cosmological 
observations at large scales. However, at second order of cosmological perturbations, the 
dipole moment entering the stress-energy tensor cannot be absorbed in an effective per- 
turbed velocity field, which means that the dipolar dark matter fluid could in principle 
be distinguished from a standard perturbed dark matter fluid. Working out the theory of 
second-order cosmological perturbations could thus yield distinctive features of the present 
model and reveal a signature of the dipolar nature of dark matter. We have particularly 
in mind effects linked with the non-gaussianity of the CMB fluctuations that are associated 
with second-order perturbations. 

F. Perturbation of the Einstein equations 

The Einstein equations at first perturbation order around the FLRW background read 

SG^ = 8tt (ST^ + J2 6T r) , ( 3 - 46 ) 

f 

where G"" = W v - \g^ v R is the Einstein tensor and where 5T" V = 8Tj£ + 5T^ is the 
perturbative part of the stress-energy tensor of the dipolar fluid given by (13.341) . The sum- 
mation runs over all the other cosmological fluids present (baryons, photons, neutrinos, . . . ) 
which are described by stress-energy tensors Tf v . Separating out the dark matter from the 
dark energy (using the link Wo = A/87r) we get 

SG^ + A SgT = 8vr (5TZ + ^ tfZf ) . (3.47) 

f 

As we have seen in the previous section, the dark matter fluid is entirely described at 
linear perturbation order by the gauge- invariant variables V, Vi and 5? (and the background 
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density p) obeying the evolution equations (I3.45P like for an ordinary pressureless fluid. 
We can thus immediately write the gauge-invariant perturbation equations in the standard 
SVT formalism (see e.g. |33|). Though these are well-known, we reproduce them here for 
completeness. For the scalar modes, we have 

A^ - 3H 2 X = 4ira 2 (p5 F + J^Pi S f) > ( 3 - 48a ) 

f 

^-$ = 87ra 2 ^p f w f a f , (3.48b) 

f 

y' + U $ = -4tt a 2 (p V + Y^ Pi (! + w d v f) » (3.48c) 

f 

H X' + (H 2 + 2H') X = 4na 2 J2p{ ( w t r f + c \ S f + f ™f Aff f) , ( 3 - 48d ) 



where we have singled out the contribution of the dipolar dark matter (cf. the variables 
V, 8p and p) from the other fluid contributions described by their background density p f , 
equation of state Wf, adiabatic sound velocity Cf, and gauge-invariant entropy perturbation 
T f . We also introduced the SVT components of the perturbative part of the anisotropic 
stress tensor, defined by 5^ = a 2 p { Wf [A ij a f + D&of + a\ j ] with A'-? = D l D j - 7 ij A/3. 
The variables Of, a\ and af are gauge-invariant because the background part of Tif vanishes 
in the case of a perfect fluid. The equations for the vector and tensor modes are 

(A + 2K) & = -16tt a 2 (p V* + J^ Pt ( l + w f) V i) > (3.49a) 

f 

® i ' + 2H® i = 8<Ka 2 ^2p i w i ai, (3.49b) 

f 

E l] " + 2HE ij/ + (2K - A)E lj = 8na 2 ^p f w f erf . (3.49c) 



f 



We highlight once more the fact that at first perturbation order, the dipolar dark matter is 
like ordinary dark matter, as can be seen from the fluid equations ( 13.451) and the Einstein 
equations (I3.48p - p.49p . Indeed, these sets of equations can be evolved without any reference 
to the dipole moment A*. 

Combining the dipolar dark matter equations (I3.45ap and f)3.45cp with the Einstein equa- 
tions (I3.48ap - fl3.48bp . we get the equation governing the growth of the dipolar dark matter 
density contrast as 

5'f + H5' F - An -pa 2 S F = 3H 2 X + 4vra 2 J^Pf ( 6 f ~ 2wf A<r f ) . (3.50) 

f 

Again, we find that the growth of structures driven by the equation (13.45cj) or equivalently 
(13.501) for the dipolar dark matter of the present model is identical with that in the standard 
CDM model at linear perturbation order. For sub-Hubble modes one can neglect the first 
term in the RHS, and we expect that the contribution of the dark matter dominates that of 
the other fluids, so we can neglect also the second term in the RHS of (13.501) . 

Interestingly, we have found in (13.291) that each of the components of the dipole moment 
obey the same equation as (13.501) but with exactly zero RHS. Recall that the dipolar dark 
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matter density contrast is defined by (13.371) as 

5 F = e F - D { \i . (3.51) 

From (13.291) we see that the internal energy due to the dipole moment satisfies the "homo- 
geneous" equation that is associated with (13.501) . viz. (recalling p = a) 

D i X , -+HD l X' i -4TTpa 2 D i X t = 0. (3.52) 

This result indicates that, in the non-linear regime, the internal energy related to the dipole 
moment may contribute significatively to the growth of perturbations (see section IIVBI for 
more comments). Finally, it is clear that the rest-mass density contrast obeys the same 
"inhomogeneous" equation, i.e. 

e F + He' F - 4vr -pa 2 e F = 3H 2 X + 4vra 2 ^p f (sf - 2w t Aa f ) . (3.53) 



IV. DIPOLAR DARK MATTER AT GALACTIC SCALES 

In this section, we shall show that, under some well motivated assumptions, the dipolar 
dark matter naturally recovers the phenomenology of MOND for a typical galaxy at low 
redshift. Such a link between a form of dipolar dark matter and MOND was the primary 
motivation of previous works [25|, [26| . We shall see that with the present improvement of the 



model with respect to [26J], thanks to the fact that the fundamental potential in the action 



now depends on the polarization field II_l = o"£j_ (instead of £j_ in the previous model |26| ) . 
the relation with MOND is straightforward and physically appealing. 

A. Non-relativistic limit of the model 

We investigate the non-relativistic (NR) limit of the dipolar fluid dynamics described by 
the equations ( 12. 21 aft and (I2.21bl) . and by the stress-energy tensor ( 12.241) . To do so, we 
consider the formal limit c — > +oo, n which is equivalent to the condition «<Cc, where v is 
the typical value of the coordinate three-velocity of the dipolar fluid. To consistently keep 
track of the order of relativistic corrections, we systematically write as O (c _n ) a typical 
neglected remainder. 

We are interested in the dynamics of dipolar dark matter and ordinary baryonic matter 
in a typical galaxy at low redshift. Let us introduce a local Cartesian coordinate system 
{ct, z 1 }, centered on this galaxy around some cosmological epoch, and which is inertial in 
the sense that it is without any rotation, nor acceleration with respect to some averaged 
cosmological matter distribution at large distances from the galaxy. Such a local coordinate 
system can be derived from the cosmological coordinate system {n.x 1 } used in section ITTH 
by posing 

ct = a(r] ) (r? - rjo) , (4.1a) 

z* = a(r, ) {x l - xl) , (4.1b) 



11 From now on, we reintroduce for convenience all factors of c and G. 
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near an event occuring at cosmological time rjo and at the galaxy's center x . In the local 
coordinate system, the metric developed at the lowest NR order reads 

97 7 

goo = -l + — + 0(c' 4 ), (4.2a) 

g 0l = O(c- 3 ), (4.2b) 

9ij =fl + ^\ 5 l3 + O (c- 4 ) , (4.2c) 

where U <C c 2 is a Newtonian-like potential. For the motion of massive (non-relativistic) 
particules we need only to include the contribution of U in the 00 metric coefficient. Thanks 
to the standard general relativistic coupling to gravity in the ij metric coefficient, the motion 
of photons agrees with the general relativistic prediction with Newtonian-like potential U. 
In the NR limit, the equation of motion (I2.21al) is readily seen to reduce to 

— - <f = -rri>v n± + o (c~ 2 ) , (4.3) 

where a 1 = dv l /dt = (d t + v^dj) v l is the standard Newtonian acceleration of a fluid in the 
Eulerian picture, v % being the coordinate three- velocity, and g l = dill the non-relativistic 
local gravitational field. Note that g % is generated by both the ordinary baryonic matter and 
the dipolar dark matter. Similarly, the equation of evolution ( I2.21bl) for the dipole moment 
reads in the NR limit [using also (14.31) ] 

^ - fll W n± = Idi (W - U ± W n± ) + fag* + O (c- 2 ) , (4.4) 

where we explicitly have d 2 ^/dt 2 = (<9 2 + a j dj + 2v j dj t + v j v k d^ k ) £]_. Notice the second 
term in the RHS which is a tidal term coming from the Riemann curvature coupling in 
(12.21bl) . Finally, the equation (12.11) reduces to the classical continuity equation 

d t a + di (av { ) = O (c- 2 ) . (4.5) 

Next, we need to be cautious about the relativistic order of magnitude of the potential 
function W appearing in the Lagrangian (12.71) . It is clear that W has the dimension either 
of a mass density or an energy density, depending of where we would reinstall the factors 
c in (j2.7p . We shall from now on assume that W is an energy density, and has a finite 
non-zero limit when c — > +oo. This will be justified when we show in ( 14.221) below that 
the coefficients W2, W3, ... in the expansion of W considered as an energy density, can be 
expressed solely in terms of G and the MOND acceleration a^ (without any c's). Therefore, 
our assumption means that we are viewing a^ as a new fundamental acceleration scale a 
priori independent from c. With such hypothesis, if we reintroduce the factors of c in 
the expression of the density r considered as a mass density and given by (I2.27al) . we get 
r = p + (W - n ± W n± )/c 2 , where p is given by (12^281) . Thus, the term (W - II ± W n± )/c 2 
becomes negligible in the formal limit c —>■ +00, and we have r = p + 0(c~ 2 ). In particular, 
we observe that the term Wo, which is linked to the cosmological constant by (restoring the 
c's and G) 

m = ^, (4.6) 
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does not enter the expression of the dipolar fluid density r, and therefore has no influence 
on the local dynamics of the dipolar dark matter in the NR limit. Our assumption that W 
has a finite non-zero limit when c — > +00 means that the cosmological constant A should 
scale with c -4 , which will be justified later when we show that A ~ Oq/c 4 . 

Thus, in the NR limit we need to consider only the mass density of the dipolar dark 
matter given by p. Now, from (12.281) we have p = a — VaII^_ which becomes when c — > +00 

p = a - diW x + O (c- 2 ) . (4.7) 

At that order the dipolar term involves only an ordinary partial space derivative. Finally, we 
get the Poisson equation in the standard way as the NR limit of the 00 and ii components 
of the Einstein equations, and find 

AU = -AttG (p h + a- <9 4 ITi) + O (c~ 2 ) , (4.8) 

where pb is the Newtonian mass density of baryonic matter. This equation can be written 
in the alternative form 

di (g i - 4tcG Hi) = -4ttG (p b + a) + O (c^ 2 ) . (4.9) 

To summarize, the equations governing the dynamics of the dipolar dark matter and the 
gravitational field in the NR limit are: the equation of motion (14.3)) . the evolution equation 
(I4.4p . the continuity equation (14.5)) and the Poisson equation (14. 9p . On the other hand, 
baryons and photons obey the geodesic equation, which means dv b /dt = d{U + (D(c~ 2 ) for 
baryons, and the standard GR formula for light deflection in a potential U for photons, 
where U is generated by (|4.8j) . 



B. The weak clustering hypothesis 

We have shown in section HTT1 that at linear perturbation order, in a cosmological context, 
the dynamics of dipolar dark matter cannot be distinguished from that of baryonic matter 
or standard dark matter. We now argue that the motion of dipolar dark matter being non- 
geodesic, its non-linear dynamics should be different. Our main motivation for the argument 
is the existence of an exact solution of the equations governing the dynamics of the dipolar 
dark matter in the NR limit. Indeed, we show in appendix |X] that, in the simple case where 
the baryonic matter is modeled by a spherically symmetric mass distribution, there is a 
solution to the equations for which the dipole moments are in equilibrium (£j_ = const), and 
at rest (v l = 0), with the internal force T % exactly balancing the gravitational field g l . In 
such a solution, the dipolar medium is uniformly distributed or more generally spherically 
symetrically distributed, and the polarization IT ± is aligned with the gravitational field g l ; 
the dipolar fluid is thus polarized. Furthermore, we show in this appendix that the latter 
solution is stable against dynamical perturbations. 

From that solution, we expect that the dipolar medium will not cluster much during the 
cosmological evolution because the internal force may balance part of the local gravitational 
field generated by an overdensity (see Fig. [2] for a picturial view of the argument). From this 
we infer that the dark matter density contrast in a typical galaxy at low redshift should be 
small, at least smaller than in the standard A-CDM scenario. Such a galaxy would therefore 
be essentially baryonic, with a typical mass density of the dipolar dark matter o rather small 
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FIG. 2: Schematic view of two worldlines of baryonic matter and dipolar dark matter. The baryonic matter 
follows a geodesic motion, u^ = 0, and therefore collapses in the regions of overdensity. Obeying the non- 
geodesic equation of motion u^ = — JF^, the dipolar dark matter is expected to have a different behavior in 
the non-linear (NL) regime. Namely, the internal force T % can balance the gravitational field g % created by 
an overdensity, in order to keep the rest mass density of dipolar dark matter close to its mean cosmological 
value, a ~ a, or at least far smaller than the baryonic one. 



compared to the baryonic one, and perhaps around its mean cosmological value a. Thus, 
the crucial hypothesis we are making (based on the solution in appendix EJ is that 



o" < Pb 



(4.10) 



or perhaps that a stays essentially at a value of the order of its mean cosmological value, 



a ~ a <C p b • 



(4.11) 



Note that for standard CDM (or baryonic matter), the density contrast between the value of 
Pcdm (or pb) in a galaxy and the mean cosmological one p cdm (or p b ) is typically of order 10 5 . 
This means that even if dipolar dark matter clustered enough so that for instance a ~ 10 3 a 
in a galaxy at low redshift, it would still satisfy the condition (14.101) . 

Note also that with this hypothesis, the non-linear growth of structures in our model will 
not be triggered by the rest mass a of dipolar dark matter (since it does not cluster much), 
but by the internal energy p- mt of the dipolar medium, which is such that p = a + p- m t and 
is explicitly given by p int = — VaII^ [recall (12.281) ]. We have seen that, at first cosmological 
perturbation order, the density contrast associated with p int reduces to —D l \ i) and obeys 
the standard evolution equation (I3.52p . We expect that at non-linear order it will take 
over the dominant role as compared to the rest mass density contrast e in the formation of 
structures. On the other hand, in the NR limit p int reduces to —diU l ± [see (14.71) ] and, as we 
shall see in the following section, will be at the origin of the MOND effect. 

We shall refer to the condition (14.101) [or even to the stronger condition (14. lip ] as the 
hypothesis of weak clustering of the dipolar dark matter fluid. Obviously, the validity of this 
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hypothesis cannot be addressed with the formalism of first-order cosmological perturbations 
in section 1111} because it is a consequence of the non-linear cosmological evolution. The 
hypothesis of weak clustering of dipolar dark matter should be validated through numerical 
N-body simulations. 

Let us thus assume that the dipolar dark matter has not clustered very much, and even 
that a might stay more or less at the cosmological mean value a (such that fidm — 0.23). 
Because of its size and typical time-scale of evolution, a galaxy is almost unaffected by the 
cosmological expansion of the Universe. Therefore, the cosmological mass density a of the 
dipolar dark matter is not only homogeneous, but also almost constant in this galaxy. Thus, 
the continuity equation f!4.5[) reduces to d- h (av l ) ~ 0. The most simple solution obviously 
corresponds to a static fluid verifying 

v*~0. (4.12) 

It is therefore natural to consider that the dipolar dark matter is almost at rest with respect 
to some averaged cosmological matter distribution. This is supported by the exact solution 
found in appendix [XJ which indicates that the dipolar dark matter in presence of an ordinary 
mass does indeed behave essentially like a static medium. Because of the internal force, the 
motion is not geodesic, and the force acts like a "rocket" to compensate the gravitational 
field and to keep the dipolar particle at rest with respect to ordinary matter (see Fig. [2]). 

C. Link with the phenomenology of MOND 

Let us now show that under the weak clustering hypothesis, the equations (I4.3l) - (l4.5p 
and (14.91) naturally reproduce the phenomenology of MOND. First of all, if (I4.12p holds, 
equation H4.3[) tells that the polarization n\ should be aligned with the local gravitational 
field g l , namely 12 

^ = fr ± W n± . (4.13) 

This proportionality relation will be the crucial ingredient for recovering MOND. 

We must now further specify the "fundamental" potential W entering the original action 
(12. 7p . In section 1111} we considered the dipolar fluid at early cosmological times, where 
the polarization field was perturbative. We shall now consider it at late cosmological times 
(around the value rjo) but still in a regime where the polarization field is weak. This will 
correspond to the outer zone of a galaxy at low redshift, where the local gravitational 
field generated by the galaxy is weak. We therefore assume that the potential W can still 
be expanded in powers of II ^ and we keep only a few terms in the expansion. Next, we 
introduce a fundamental acceleration scale a to be later identified with the MOND constant 
acceleration whose commonly accepted value is a$ — 1.2 x 10~ 10 m/s 2 [8(. Associated with 
Qq we can define a fundamental surface density scale 



2ttG 



(4.14) 



whose numerical value is X ~ 0.3 kg/m 2 ~ 130 M /pc 2 . The numerical value of S is 
close to the observed upper limit of the surface brightness of spiral galaxies — the so-called 



12 From now on, we no longer indicate the neglected remainder terms O(o 2 ). Furthermore we assume for 
the discussion that (|4. 12[) is exactly verified, i.e. v l = 0. 



2(3 



Freeman's law which is seen as an empirical evidence for MOND [8(. We now assert that the 
expansion of W when n_i_ — > is physically valid when the condition II j_ <C S is satisfied. 
As will become obvious, this condition can equivalently be written g <C a , where g = \g l \ is 
the norm of the local gravitational field of the galaxy, and this will correspond to the deep 
MOND regime (see Fig. ED - With respect to the expansion (13. 14|) already used in cosmology, 
we shall be able to add an extra term. We now write this expansion, for II_l <C E, as 



w(n ± ) 



w + ~ m ni + ~ m n\ + o 



m 



(4.15) 



iw 2 ni + i; 

The first term Wo is related to the cosmological constant A through ( 14.61) . We now show 
that the next two coefficients W2 and W3 are uniquely determined if we want to recover the 
phenomenology of MOND. Indeed, by inserting (14.151) into the relation (14.131) we obtain 
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m + - w 3 n_L + o 



m 



(4.16) 



which can be inverted to yield the polarization as an expansion in powers of (the norm of) 
the gravitational field. Anticipating that W2 S ~ ao, this expansion will be valid whenever 
g <^ do. We obtain 

m 



if, 



1- - 



2 W ; 9 + 



(4.17) 



Next, following the conventions of [25j, |26|, we introduce the coefficient of "gravitational 
susceptibility" \ of the dipolar medium through 



IF, 



X 



AttG 



Inserting that definition 13 into the FHS of the Poisson equation (14.91) . we find 

d i [(i+ x )g i ] = -47TG (p b + a) . 



(4.18) 



(4.19) 



Finally, invoking our hypothesis of weak clustering (14. 10J) . or (14. lip in the more extreme 
variant, we can neglect the mass density a of the dipole moments with respect to the 
baryonic one, so we obtain the MOND equation which is generated solely by the distribution 
of baryonic matter as 39J 

d i (ng i ) = -4irG Ph . (4.20) 

The MOND function \x is related to the susceptibility coefficient by fx — 1+vand can actually 
be interpreted as the "digravitational" coefficient of the dipolar medium [25j. Again, let us 
stress that in this model we do have some distribution of dark matter a in an ordinary sense, 
but we expect its contribution to become negligible in galactic halos at low redshifts (after 
cosmological evolution), so that the MOND fit of rotation curves of galaxies is unaffected 
by this "monopolar" dark matter. 14 The MOND effect is due to the dipolar part of the dark 
matter given by the internal energy p mt = —^11^. 



13 Note that this definition is valid in both MOND and Newtonian regimes whenever the polarization is 
aligned with the gravitational field. 

14 However, at the larger scale of clusters of galaxies the monopolar part of the dipolar medium a may play 
a role to explain the missing dark matter in MOND estimates of the dynamical mass [8|, |40( • Note that in 
the present model, the motion of photons, needed to interpret weak-lensing experiments, is given by the 
standard general relativistic prediction; see (|4.2p with potential U solution of the MOND equation (|4.20p . 
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Now, from astronomical observations we know that the gravitational susceptibility x i n 
the deep MOND regime g <C a® should behave like 

X = -l + - + o(U] 2 ). (4.21) 



«o 



«n 



The fact that x should be negative was interpreted in the quasi-Newtonian model [25 
as an evidence for gravitational polarization -- the gravitational analogue of the electric 
polarization in dielectric media. By inserting (I4.2ip into (I4.18p . and comparing with the 
prediction of our model as given by ( 14.1 7ft . we uniquely fix the unknown coefficients therein 
as 

W 2 = 4vrG , (4.22a) 

(~i2 

W 3 = 32vr 2 — . (4.22b) 

a 

This, together with Wo fixed by (14.61) . determines the potential function W up to third order 
from astronomical observations. As we see, the MOND acceleration ao enters at third order 
in the expansion, and therefore does not show up in the linear cosmological perturbations 
of section IHII At third order, the potential W deviates from a purely harmonic potential, 
and ao can be seen as a measure of its anharmonicity. 

To express W in the best way, we prefer using the surface density scale £ = a /2iiG 
rather than the acceleration scale ao- To do so, we must introduce a purely numerical 
dimensionless coefficient a to express the cosmological constant A (which is positive and has 
the dimension of an inverse length squared) in units of ao/c 4 , and we pose 
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A = 3a^-^J . (4.23) 

The definition of a is such that a^ = aa® represents the natural acceleration scale associated 
with the cosmological constant, and is already given by (11.11) as a a = a/A/3c 2 /27t. Then, 
the cosmological term (14.61) becomes Wo = 6n 3 CE 2 a 2 , and we obtain 



W = 6ttGS 




E 



(4.24) 



In the present model there is nothing which can give the relation between A and ao, hence 
a is not determined. However, if the dipolar fluid action (12. 7p is intended to describe at 
some macroscopic level a more fundamental theory (presumably a QFT), we expect that 
the potential W should depend only on certain more or less fundamental constants, and 
some dimensionless variables built from "fundamental fields" . Introducing the dimensionless 
quantity x = IIj_/E, we can rewrite (14.241) as W = QtiG £ 2 w(x), where 

w ( x ) = a 2 7T 2 + -x 2 + -x 3 + C(x 4 ) (4.25) 

3 9 

represents some "universal" function coming from some fundamental albeit unknown physics. 
Therefore, we expect that the numerical coefficients in the expansion of w(x) should be of 
the order of one or, say, 10. In particular, it is natural to expect that a should be of the 
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order of one (to within a factor 10 say), and we deduce from (14.231) that the magnitude of A 
should scale approximately with the square of the MOND acceleration, namely A ~ a^/c 4 . 



The numerical coincidence between the measured values of A and a® is well-known 16 
The observed value of the cosmological constant is around A ~ 0.12 Gpc -2 [33J which, 
together with a^ ~ 1.2 x 10~ 10 m/s 2 , corresponds to a value for a which is very close to one: 
a ~ 0.8. Thus a is very close to the scale oa associated with the cosmological constant, 
which is related to the Gibbons-Hawking temperature Tqh = ha^/kc derived from semi- 



classical theory on de Sitter space-time [411 ] . From the previous discussion, we see that the 
"cosmic" coincidence between A and a^ has a natural explanation if dark matter is made of 
a medium of dipole moments. 

D. The Newtonian regime 

For the moment, we looked at the explicit expression of the potential function W in the 
MOND regime g <^i clq. We would also like to get some information about this function 
in the Newtonian regime g ^> a®. To do so, we first derive the general expression of the 
gravitational susceptibitity coefficient x- Here we assume that the MOND function n = 1 + x 
is always less than 1. This implies x < an d thus using (14.131) and (I4.18P we must have 
W n > (where we recall that W n = dW/dn_|_). Taking the norm of (I4.13P we get 
g = W n (n_i_). Next, we introduce the function ©(g) which is the inverse of W n (n_i_), i.e. 
satisfies 

g = W n± (U ± ) <{=► IL ± = e(g). (4.26) 

According to (I4.18p . the susceptibility \ * s then given as the following fonction of the grav- 
itational field g, 

x(g) = -4nG^. (4.27) 

9 

This is the general relation linking x (or equivalently the MOND function fi = 1 + x) to the 
potential function W in the dipolar action (12.71) . Of course, in the present model W is to 
be considered as more fundamental than x which is a derived quantity. 

In the Newtonian regime g 3> ao, the MOND function fi should tend to one, so that 
X vanishes in this regime. To discuss more concretely this condition, we assume that in 
the formal limit g — > +oo, the gravitational susceptibility behaves as x ~ g~ J , with 7 a 
strictly positive real number. More precisely, it should behave like x ~ ~ e (5 , / a o) _7 5 where 
e is a strictly positive real number. Beware that even if this power-law behavior is a simple 
assumption, nothing garanties that it is verified. Then, when g — > +00, we get from (I4.26P 
and (flbrp that 

n^-A^ 1 " 7 , (4.28a) 

W~^—^Ag 2 -~< + K, (4.28b) 

2-7 

where A = eal/^itG > and k is an integration constant. We have to distinguish several 
cases, depending on the value of the exponent 7: 

(i) If < 7 < 1, then both the polarization n^ and the potential W diverge. This would 
corresponds to the curve (a) of Fig. HI 

29 




Newtonian regime 
(a) 



n, 




n, 



FIG. 3: The minimum of the potential function 
VV(II_i_), reached when II j_ = 0, is a cosmological 
constant A. Small deviations around the mini- 
mum, corresponding to II j. <C S = a /2nG, de- 
scribe the MOND regime g <C an- 



FIG. 4: The potential W as a function of the po- 
larization II j_ for different asymptotic behaviors 
of the gravitational susceptibility \ m the New- 
tonian regime </ > aq. The arrows indicate the 
direction of increasing gravitational field g. 



If 7 = 1, the polarization n_i_ tends to a maximum "saturation" value IT r 
the potential W equals the constant k. See curve (b) in Fig. HI 
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If 1 < 7 < 2, the polarization goes to zero while the potential diverges to — oo like 
a power law. This implies that W cannot be a univalued function of IIj_. Therefore, 
there must exist two branches corresponding to the Newtonian and MOND regimes. 



If 7 = 2, according to (14.28b|) the potential diverges to — oo logarithmically, i.e. W ~ 
—A lng, while the polarization still vanishes. Same conclusions as in case (iii) apply. 

(v) Finally, if 7 > 2, the polarization goes to zero while the potential tends to k. Same 
conclusions as in (iii) apply. 

If we believe that the potential W represents a fundamental function in the action, and that 
our model should strictly speaking be valid in a Newtonian regime (and not being merely 
valid in the MOND regime), we should a priori expect that W is a univalued function of n^. 
Then, the susceptibility coefficient should be like \ ~ g~ 7 with < 7 ^ 1 in the Newtonian 
regime. This would mean that the MOND function /1 behaves like 



/' 
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(4.29) 



with < 7 =Sj 1. Such rather slow transition of \x toward the Newtonian regime is consistent 
with the recent results of |42j who fitted the rotation curves of the Milky Way and galaxy 



NGC 3198, and of 43] who fitted 17 early-type disc galaxies, and concluded that the New- 
tonian regime is rather slowly reached. For instance, the authors of [42j, |43|, |44| agreed that 
7 = 1 yields a better fit to the data than 7 = 2. 
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The case 7 = 1 (curve (b) in Fig. [4]) corresponds to an interesting physical situation in 
which the dipolar medium saturates when g — > +00, at the maximum value Il max = A, or 

n max = ^S, (4.30) 

where S is the surface density scale ( 14.141) . In this saturation case, the gravitational suscep- 
tibility coefficient behaves as 

X-~e^. (4.31) 

9 

However, let us remind that such a slow transition from MOND toward the Newtonian 

regime is a priori ruled out by Solar System observations. Indeed, according to the MOND 

equation, a planet orbiting the Sun feels a gravitational field g obeying (l + x)9 = <7n, where 

#n is the Newtonian gravitational field. Hence, if \ scales like g^ 1 when g ^> a like in 

f )4.3ip . the gravitational field experienced by planets will involve a constant supplementary 

acceleration directed toward the Sun (i.e. a "Pioneer-type" anomaly) given by 

g~g^ + ea . (4.32) 

Of course it is striking that the order of magnitude of the Pioneer anomaly is the same as 
the MOND acceleration clq. Unfortunately, the presence of a constant acceleration such as 
in (14.321) should be detected in the motion of planets, and this is incompatible with current 
measurements (see e.g. 45|, |46j for a discussion). 



Despite the fact that a slow transition to the Newtonian regime (like for example the 
case 7 = 1) seems to be favored by observations at the galactic scale [42j, |43|, |44j, it does 



not seem to be viable when extrapolated up to the scale of the Solar System. In our model, 
we found that such a behavior is the result of our belief that the "fundamental" function 
W be univalued. In this respect, the validity of the model should be limited to large scales, 
from the galactic scale up to cosmological scales, i.e. in a regime of weak gravity. At smaller 
scales the description in terms of a single univalued function W should break down. But 
our model being an effective one, or even a phenomenological one, the question of whether 
the potential W is univalued or not remains an open issue. 

V. SUMMARY AND CONCLUSION 

In this paper, we proposed a model of dark matter and dark energy based on the concept 
of gravitational polarization of a medium of dipole moments. The dynamics of the dipolar 
fluid is governed by the Lagrangian (12.71) in standard general relativity, and constitutes 
a generalization of the previous model [26j. Namely, this Lagrangian involves a potential 
function W, describing at some effective level a non-gravitational internal force influencing 
the dynamics of the dipolar fluid, and which depends on the polarization field or density of 



dipole moments n^ = o~£,± instead of merely the dipole moment itself £j_ in the model [26 
This new form of the potential permits recovering in a most elegant way the phenomenology 
of MOND in a typical galaxy at low redshift. In addition, we show that the model naturally 
contains a cosmological constant A. 

We proved in section IIHI that whithin the framework of the theory of first-order cos- 
mological perturbations, the dipolar fluid behaves exactly as standard dark energy (i.e. a 
cosmological constant) plus standard dark matter (i.e. a pressureless perfect fluid). Thus, 
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our model is consistent with the cosmological observations at large scales. In particular, 
it leads to the same predictions as the standard A-CDM model for the CMB fluctuations. 
However, at second order in the cosmological perturbations, we expect that the dipolar 
dark matter should differ from a perfect fluid because of the influence of the internal force 
resulting in a non-geodesic motion. The model could thus be checked by working out the 
second-order cosmological perturbations and comparing with CMB fluctuations (notably the 
effects linked with the non-gaussianity) . 

The dynamics of the dipolar dark matter being different from that of standard dark 
matter (at the level of non-linear perturbations), we expect the "monopolar" part of the 
dipolar dark matter not to cluster much during the cosmological evolution. We call this 
expectation the hypothesis of "weak clustering" . It is supported by an exact solution worked 
out in appendix |X] for the dynamics of dipolar dark matter in the non-relativistic limit and 
in spherical symmetry. In this solution, the internal force balances the local gravitational 
field produced by a spherical mass, so that the dark matter remains at rest with respect to 
the central mass. The weak clustering hypothesis should be checked via N-body numerical 
simulations. Under that hypothesis, we show that the Poisson equation for the gravitational 
field generated by the baryonic and dipolar dark matter reduces to the MOND equation in 
the regime of weak gravitational fields g <C ao- Our model of dipolar dark matter therefore 
naturally explains all the successes of the MOND phenomenology. 

To achieve this result (in section HVl) we have to adjust the fundamental potential W in the 
action. We find that it should be given by an anharmonic potential, the minimum of which, 
reached when U± = 0, being directly related to the cosmological constant A. It is tempting 
to interpret A as a "vacuum polarization" of some hypothetical quantum field, when the 
"classical" part of the polarization U± — > 0. The expansion around that minimum is fine- 
tuned in order to recover MOND. In particular, we show that the MOND acceleration ao 
parametrizes the coefficient of the third-order deviation of W from the minimum. Although 
fine-tuned to fit with observations, this potential function W offers a nice unification between 
the dark energy in the form of A and the dark matter in the form of MOND (see Fig. [3]). 
A consequence of such unification is that the cosmological constant should scale with the 
MOND acceleration according to A ~ ag/c 4 . This scaling relation is in good agreement with 
observations and has a very natural explanation in our model. 

To conclude, we proposed to modify the matter sector rather than the gravity sector 



as in modified gravity theories [14|, [18|, [23|, [24|]. Namely, we investigated a model of dark 



matter, but of such an exotic form that it naturally explains the phenomenology of MOND at 
galactic scales. Furthermore, that form of dark matter has a simple physical interpretation 
in terms of the well-known mechanism of polarization by an exterior field. More work is 
necessary to test the model, either by studying second-order perturbations in cosmology, 
or by computing numerically the non-linear growth of perturbations and comparing with 
large-scale structures. 
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APPENDIX A: DARK MATTER IN A CENTRAL GRAVITATIONAL FIELD 

We investigate the dynamics of the dipolar dark matter fluid in presence of a spherically 
symmetric mass distribution of ordinary baryonic matter in the NR limit c — ► +00. The 
equations to solve are the equation of motion (14. 3p . the equation of evolution (14.41) . the 
continuity equation (14 .5p and the Poisson equation for the gravitational field (14.91) . Let us 
rewrite those equations here for convenience: 15 

-^=g-J=, (Ala) 

d t a = -V • (av) , (Alb) 

V • g = -AttG (a + p b - V • n) , (Ale) 

|^f = T + ^ V (W - n W) + (£ • V) g , (Aid) 

where the internal force reads J- = II W, with II = 11/11. 

Our aim is to solve the equations (lAip in the special case where the baryonic matter is 
modeled by a time- independent spherically symmetric distribution of mass Ph{f), say with 
compact support. Let us show that there is a simple solution to such a set of equations, in 
the case where 

v = , (A2a) 

a = a (r) , (A2b) 

which corresponds to a static fluid whose mass distribution is time-independent and spher- 
ically symmetric. We denote such particular solution with a lower index 0. From (1A2J) we 
observe that the continuity equation (lAlbl) is immediately satisfied. In such a solution, 
according to (lAlap the internal force balances exactly the gravitational field, i.e. J-q = gr 
(this is somewhat similar to the case of a non-rotating star in hydrostatic equilibrium, where 
the pressure gradient plays the role of the internal force). We deduce that the polarization 
field n = cro^o is aligned with the gravitational field go- Hence, from equation (lAlcft both 
n and go are radial. We shall pose go = —go(r,t) e r and n = — n (r, t) e r , where in our 
notation g > and n > 0. 

Furthermore, let us show that in addition the polarization field is practically in "equilib- 
rium", i.e. n is independent on time t, and so is go. We replace go by the explicit expression 
of the internal force To = FToWq into the evolution equation (lAldp . use (!A2aP and get 

djiio - aoW'o rio = v (Wo - n Wo) + (n • V) (no W' Q ) . (A3) 

Here IIo = n /n = — e r , and we introduced the shorthand notation Wq = W"(n ). Now, 
it turns out that the RHS of this equation vanishes in the special case where the polarization 
field is radial, hence we get 

«9 t 2 n = aoW'o . (A4) 



15 In this appendix, we adopt 3-dimensional notations with boldface vectors, e.g. T = {T l ). We also remove 
the subscript _L from the variables £±_ and II j_ for notational simplicity. The derivatives of the potential 
W with respect to its argument II will be denoted with a prime, e.g. W' = W n = dW/dll. 
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In order to determine the time evolution of Ho, an explicit expression for the potential W is 
in principle required. However, we saw in section IIVCI that the potential W only depends 
on the polarization II and the constants a® and G. The only time-scale one can build with 
a , G and a is the dipolar dark matter self-gravitating time-scale r g = (Tr/Gao) 1 / 2 , or 
equivalently, in terms of frequency, uo 2 = AnGaQ. Therefore, the polarization n can only 
evolve on this time-scale. For instance, in the MOND regime g <C ao, we have at leading 
order W = 4iiG n , hence ( 1A4D reduces to 

d 2 n = ^ 2 n . (A5) 

The most general solution of this equation is a linear combinaison of hyperbolic coshu; g £ and 
sinhtUgt. For a "monopolar" dark matter mass density <7o of, say, the mean cosmological 
value a ~ 10 -26 kg/m 3 [in agreement with our weak clustering hypothesis (14. lip ], the typical 
time-scale of evolution of n will be larger than 6 x 10 10 years. This is large enough to neglect 
any time variation of il with respect to a typical orbital time-scale in a galaxy. Our solution 
is therefore given by 

n = -Ilo(r) e r , (A6) 

together with (1A2I) . The dipole moments are at rest and in equilibrium. The explicit function 
Ilo(r) is determined from the radial gravitational field go (r) as 16 

Mr) = © Mr)) , (A7) 

where O(g'o) denotes the inverse inverse function of W"(n ) following the notation (14.261) . 
The gravitational field go(r) is determined by the Poisson equation (I Aid) as 



90 -4,Gn =^|M, (A8) 

where M (r) = 4n J Q dss 2 [pb(s) + o"o(s)] is the mass enclosed within radius r. 

The existence of this physically simple solution represents a notable progress compared 
to the more complicated solution found in the previous model [26] (see section IV there). 
Such a solution is quite interesting for the present model because it indicates that during 
the cosmological evolution (at non-linear perturbation order) the dipolar dark matter may 
not cluster very much toward regions of overdensity. Most of the effect will be in the dipole 
moment vectors which acquire a spatial distribution. This is our motivation for the "weak 
clustering" assumption (I4.10I) - (I4.11I) stating that a <C p b , which was used in section IIVCI 
to obtain MOND. In the present case, neglecting er with respect to pb i n the RHS of 



we recover the usual MOND equation generated by the baryonic density only. This being 
said, such an appealing solution may be physically irrelevant if the spherically symmetric 
configuration appears to be unstable with respect to linear perturbations. This motivates 
the following study of the stability of the previous solution. 

Consider a general perturbation of the background solution, namely 

a = a + 5a , (A9a) 



16 Note that if in this solution the polarization field ilo(r) = <Jo(r)£,o(r) is determined, the density ao(r) 
and dipole moment £o(0 are not specified separately. For instance, the density could be at the uniform 
cosmological value a so that £o( r ) = IIo(r)/o r . This degeneracy of <Jo(r) is an artifact of our assumptions 
of spherical symmetry and staticity. 
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V 

n 



5v , 

n + su . 



(A9b) 
(A9c) 



We have also g = g + 5g and J- = J- + 5^F, where the expression of the perturbed force 
in terms of the perturbed polarization explicitly reads 



sj= = vitf (n • sit) n + w 






n .-)n 



(A10) 



Assuming a Fourier decomposition for any perturbative quantity SX, we write for a given 
mode of frequency u and wave number k, 



SX(x,t) = 5X(k,cu)e 



i(k-x-ujt) 



(All) 



We want to find the relation between k ■ e r and u, the so-called dispersion relation, which 
contains all the physical information about the behavior of the generic perturbation (lAllft . 
Introducing this ansatz into (lAip . and simplifying the resulting equations by making use of 
the background solution, we find 



5v 
5a 

Sg 



(Sg - 8J=) , 
<T k ■ Sv — i 5v ■ Vo" ) 
4ttG *— (6cr-ik- Sit) . 



1 



(Al2a) 
(A12b) 
(A12c) 



These algebraic expressions can be combined to express 5a, Sg and Sv in terms of 5U only. 
After some algebra, we get from the evolution equation ( lAldl) a relation expressing the 
perturbed polarization field 511 = ao5£ + 5a £o as 



u 2 5U 



5a 



lO) 



"\ 



uj- — n + — (Sv ■ Vct ) n - iu (Sv ■ v) n + (n • sit) v (n w ' 
+ n w ' (n • sit) i fe - (i fe - n ) 5 & - (<m • v) s- - ^o ^ . 



(A13) 



When replacing 5cr, 5gr, 5v and 5jF into (IA13[) we obtain a master equation for the perturbed 
polarization 5H which is quite complicated. Given the complexity of the problem, we restrict 
our analysis to the simplest modes in a spherically symmetric background, namely those 
propagating radially. We shall thus write k = k e r , and study successively the transverse 
and longitudinal perturbations. 

Firstly, let us consider a transverse perturbation SH, i.e. one which satisfies 511 ■ e r = 0. 
Projecting the master equation (1A13J) in the direction of 5U, we get that 



uj 2 + W' n [^-- 
£o r 



5U = 



(A14) 



which simply states that no transverse perturbations propagating radially are allowed, i.e. 
5U = 0. Consider now the case of a longitudinal perturbation 5U = —5U(r, t) e r , where 5U 
can be positive or negative (with our convention the norm of II reads II = n + 5U) , and 
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represents the arbitrary amplitude of the applied linear perturbation. After some lengthy 
calculations, we get the dispersion relation 



,2 



O"0 V w g 



1 + 



(AtiG - Hff) <9 r n 
LU 2 + a W[; + Tl d r W[; 



(A15) 



Notice first that, as the wave number k is purely imaginary, such a perturbation cannot 
propagate. Secondly, the stability of the background solution with respect to this perturba- 
tion is related to the sign of k/i, so an explicit expression for the potential W is required to 
conclude. Such an expression is available in the MOND regime go <C a® using the expansion 
(I4.24p . Assuming the MOND equation with a (baryonic) point mass M for simplicity, i.e. 
equation (1A8J) with pb = M S(x) and negligible <To, we find that the dispersion relation can 
be recast at the leading order in the form 

a ^ + 2u; g 2 u; 2 H-a;|(^-2c4)' l J 

where io\ = GM/r 3 denotes the Keplerian frequency. We now turn to the analysis of the 
two factors in (1A16I) . namely the ^-dependent and o"o-dependent ones. 

At a given distance r from the center of the galaxy, the tu-dependent factor becomes very 
large in the vicinity of the resonant frequency 

V2u K -u g ). (A17) 



u R — u g 

But we are restricting our attention to perturbations in the MOND regime where go "C Oo, 
which means at distances r from the galactic center that are far larger than the MOND radius 
ru = a/GM/oo, or equivalently at Keplerian frequencies wk "C Wm with u^ = GM/r^. For 
a typical galaxy of mass M ~ 10 11 M , and a "monopolar" dark matter mass density 
around the mean cosmological value, i.e. o"o ~ o 7 — 1CF 26 kg/m 3 , we find by reporting the 
constraint cu K <^ cu M into (1A17J) the upper-bound u>^ <C -y/^tUgCUM, which gives numerically 
cur <^ 10 -17 s _1 . Because perturbations with a typical time scale 2ti/u> >2x 10 10 years 
are out of the present scope, the u;-dependent part of ( 1A16I) reduces to a numerically small 
factor. 

Finally, we consider the o"o-dependent part of (1A16|) . Consistent with the "weak clustering 
hypothesis" f l4.10l) - fl4.11l) . we are expecting the background density profile <7o to be almost 
homogeneous. Thus, the factor d r ao/cro will be of the order of the inverse of the characteristic 
length scale L of variation of <To assumed to be far larger than the typical size £ of the galaxy, 
which implies \k ■ x\ < tjL ~ in flAllj) . A longitudinal perturbation would therefore keep 



oscillating at the frequency uj without propagating, and we conclude that it would be stable. 
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